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Abstract 

We consider the Ising model on a cylindrical lattice of L columns, with 
fixed-spin boundary conditions on the top and bottom rows. The sponta- 
neous magnetization can be written in terms of partition functions on this 
lattice. We show how we can use the Clifford algebra of Kaufman to write 
these partition functions in terms of L by L determinants, and then further 
reduce them to m by m determinants, where m is approximately L/2. In this 
form the results can be compared with those of the Ising case of the superinte- 
grable chiral Potts model. They point to a way of calculating the spontaneous 
magnetization of that more general model algebraically. 
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1 Introduction 

Onsager calculated the partition function of the two-dimensional square-lattice 
Ising model in 1944. p] He did so by observing that the row-to-row transfer 
matrices generated a finite Lie algebra. This method was refined by Kaufman 
in 1949 [2], who showed that the matrices could be expressed in terms of a 
Clifford algebra. By taking the thermodynamic limit of a large lattice, one 
obtains the free energy per site. 

The calculation of the spontaneous magnetization is much harder than that 
of the free energy. Onsager first announced his famous result for the spon- 
taneous magnetization at a conference at Cornell University in August 1948 
and later at a IUPAP conference in Florence. [3j d] Onsager never published 
his derivation - the puzzle was finally solved by Yang in 1952. [5 J 

Here we wish to re- visit this problem, with a view to seeing if the techniques 
can be generalized to the TV-state superintegrable case of the solvable chiral 
Potts model, which has properties very similar to the Ising model and reduces 
to it when N = 2. The author has previously used large-lattice functional 
relations and analytic methods to obtain the spontaneous magnetization of 
the general solvable chiral Potts model[6l[7j, but it would still be interesting 
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to have an algebraic method that may give more insight into the properties 
of the model on a finite lattice. 

In section 2 we define the Ising model on a cylindrical lattice of L columns, 
with fixed-spin boundary conditions on the top and bottom rows. We show 
how the magnetization can be expressed, firstly as the ratio of partition func- 
tions, and consequently as a matrix element between the two ground-state 
eigenvectors of the transfer matrix T. In section 3 we introduce the hamilto- 
nian Ti which commutes with T and define "hamiltonian partition functions" 
Z, W. These are limiting cases of the usual partition functions and are slightly 
simpler to work with. This Z corresponds to the usual partition function, W 
to the partition function with an extra factor a±, where a\ is some spin deep 
inside the lattice. 

In section 4 we set up the apparatus of the Clifford algebra and use it 
to write Z, W as square roots of L by L determinants. Up to this point 
our calculation parallels that of Yang. [5 J Montroll et a/[3] also calculated the 
spontaneous magnetization using the more combinatorial pfaffian method and 
Szego's theorem on Toeplitz matrices. [8] In section 5 we marry Yang and Mon- 
troll et aVs techniques by using Szego's theorem to evaluate the appropriate 
large- L limit of W and hence derive Onsager's famous formula for the spon- 
taneous magnetization. 

In section 6 we show how the expressions for Z, W can be reduced to linear 
expressions in determinants of dimension \{L — Q)/2], where Q = or 1. Then 
Z is precisely the hamiltonian limit of the corresponding partition function 
of the two-state sup erintegr able chiral Potts model. In section 7 we write the 
result for W in terms of two orthogonal matrices B + , B_ . 

In section 8 we comment on whether our result for W can be generalized 
to the iV-state superintegrable chiral Potts model. 



The model is defined on the square lattice C, rotated through 45°, with M + 1 
horizontal rows, each containing L spins, as in Fig. 1. 

We impose cylindrical boundary conditions, so that the last column L is 
followed by the first column 1. At each site i there is a spin Cj, taking the 
values +1 or —1. The spins in the bottom row are fixed to have value a = dbl, 
those in the top row to have value +1. Adjacent spins (Ji,(Jj on southwest to 
northeast edges (with % below j) interact with Boltzmann weight ex.p(Kai(Tj); 
those on southeast to northwest edges with weight exp(Kaiaj). 

An important associated parameter is 



The system is ferromagnetically ordered if K, K are positive real and 



To align with our notation for the chiral Potts model[6j[7] we are defining the 
RHS of (12. lj) to be k', rather than k as in Onsager.pJ eq. 2.1a] 



2 The model 




(2.1) 



< k' < 1 



(2.2) 



2 




Figure 1: The square lattice C turned through 45°. 



Partition function 

The partition function, which depends on a, is 

Za = II eyLp(KaiCrj) J| exp(K<7iCTj) , (2.3) 

°~ (id) (i,j) 

the products being over all edges of the two types. The sum is over all values 
of all the free spins. The partition function can be written as 

Z a = u\T M u + , (2.4) 

where T is the row-to-row transfer matrix, with elements 

L 

T a y = HexpiKa.a'^expiKai^) , (2.5) 
i=i 

a being the set of all spins o\, . . . , ax in one row, and a' being the set in the 
row above. Thus T is a 2 L by 2 L matrix. The vector u a is of dimension 2 L , 
with entries 

(u a ) a = 1 if ax = ■ ■ ■ = (Tl = a , 

= otherwise . (2-6) 

Two transfer matrices T,T', with different values of K and K, commute 
provided they have the same value of fe'.[9l section 7.5] 

Here denoted the transposed conjugate of u. However, all our matrices 
are real, and T commutes with its transpose, so all the eigenvectors we shall 
discuss are also real: the complex conjugation is unnecessary. 

We shall also need the matrices Si, ... , Sl, Ci, . . . , Cl, defined by 

L 

(Sj)a,a> = <Tj 5(a n ,a' n ) , (2.7) 
n=l 
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(CjW = 6(^,-^)11* 5(a n ,a' n ) , (2.8) 

n=l 

where the * means that the term n = j is excluded from the product, If /, S, C 
are the two-by-two matrices 

'-(;;) C-°o •—(:;) ■ (2 - 9) 

then we can write Sj , Cj as 

the S, C on the RHS being in position j in the sequence of L factors. 
Spontaneous magnetization 

Take a = +1, so all top and bottom boundary spins are +. Let be a site 
deep within the lattice. The expectation value of its spin ctq is 

M = (cto) = ^^^oII^P^^OlI^P^^) • ( 2 - 10 ) 

a (iJ) 

We take the limit when the lattice is infinitely large, so L, M — > oo, and site 
is infinitely far from the boundaries. 

The products are unchanged by negating all spins cjj, so if we imposed 
toroidal boundary conditions, then it would be true that 

(ffo) = -Wo) (2.11) 

and this would imply that Ai = 0. At high temperatures (k' > 1), this 
is true also for our fixed-spin boundary conditions when we take the large- 
lattice limit. However, at lower temperaturers (0 < k' < 1) the system has 
ferromagnetic long-range order and "remembers" the boundary conditions 
even in the limit of deep inside a large lattice, so that 

< M < 1 . (2.12) 

If is in column i of row j + 1, then we can write ([2.100 as 

M = u ] + T j SiT M - j u + /u ] + T M u + (2.13) 

Because of our cylindrical boundary conditions, we can cyclically permute the 
columns of £, so the RHS must be independent of i. 

The sub-spaces 

The operator that negates all the spins in a row is 

R = dC 2 ---C L . (2.14) 



4 



We can divide the full 2 L -dimensional space into two orthogonal subspaces 
V+, V_ such that 

Rv = rv if v £ V r , r = + or — (2-15) 

The operator R commutes with T and Ti., so its eigenvectors lie either in V+ 
or V_. Let A a be the largest eigenvalue of T in V a , i>a the corresponding 
eigenvector normalized so that ipa^a = lj and define 



v+ = (u + + u_)/V2 , v- = (u + -U-)/V2 . (2.16) 
Then T 3 v a £ V a for all j and 

vlT M v b = if a + b , (2.17) 
vlT j SiT M -iv b = if a = b , (2.18) 

so (|2.13p becomes 

M = -± ^ = = . 2.19 

vlT M v+ + v ] T m v 



Asymptotic degeneracy 

If j,M — j are large we can take T J = A^ + ifj + ^ + or A~Lip_il)_ (depending on 
the sub-space in which T is acting). Similarly for T M ~i and T . Thus the 
four matrix elements in (|2.19p will be proportional to AiA_ A jl _A^ r_ ' ? , 
A^f, A_ , respectively. However, for < < 1, A + and A_ are asymptotically 
degenerate, in that for L large 

A+/A_ = 1 + 0(e- £ <) , (^+)/(^-) = l + 0(e~ iC ) , (2.20) 

choosing the signs of the ip a appropriately. Here Q is independent of L and is 
a measure of the interfacial tension between the two phases. section 7.10] 

That this should be so can be seen from a low-temperature series expansion 
of the eigenvalues A and eigenvectors ifi of T (or more conveniently the 7i of 
the next sub-section) in increasing powers of k'. 

One can start from a single zero-temperature configuration in which either 
all spins are +1, or all are —1 («+ or u_). The expansion then proceeds 
identically for each choice (going from one choice to the other by negating 
all spins) until one reaches terms of order k' L , when one has to consider the 
opposite state (all spins —1 or all +1) from the original. This has the same 
eigenvalue for k' = 0, so naive eigenvalue perturbation theory fails. Then and 
only then does one have to decide whether to symmetrize or anti-symmetrize 
the eigenvector with respect to R. 

If are the corresponding two near-eigenvectors, then up to this 

order the eigenvalues are the same and 

ul<j) + = v} + (j>_ = . (2.21) 

The actual eigenvectors of T are obtained by symmetrizing: 

V>+ = (</»+ + <U/V2 , ^- = (<P+ - 0-)/V2 • (2.22) 
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The relations ([gjgUJ) then follow, using (fXTU|) . 

Hence in the limit of j, M — j,L all large we can take A + = A_, and ()2. 19|) 
becomes 

M = il>\Sii>_ • (2-23) 



3 The associated hamiltonian H 

Consider the limit when K becomes small and K becomes large while k' 
remains fixed. Then to first order in K we obtain from (|2.5p 

e" L7? T = 1 - KH = e~ KH , (3.1) 
where I is the 2 i -dimensional identity matrix, 

n = Ho + k'Hi , (3.2) 

where 

L L 
i=l i=l 

and Sl+i = S\. It follows that this hamiltonian 7i commutes with the transfer 
matrix T and with R. Its ground state eigenvector (the one corresponding to 
the most negative eigenvalue of H , and to the largest eigenvalue of T) in the 
subspace V a will therefore be tp a - 

When k 1 = , the hamiltonian reduces to the diagonal matrix TCq and 
we see that the ground state eigenvectors are indeed u + and u_, both with 
eigenvalue —L. it is convenient to define the closely related matrix 

L 

J = H + LI = J2{Z-$iSi+i) ■ (3-4) 
i=l 

This has minimum eigenvalue 0, the corresponding eigenvectors being v+ and 
v_. 

Hamiltonian partition functions 

We could continue to look at the four matrix elements in (|2.19p , or attempt to 
evaluate directly the result (I2,23h , We prefer to follow an intermediate path 
and to consider the expressions 

Z+(a) = vle- an v + , Z_(P) = vlz~ m v_ , (3.5) 

W(a,p,x) = vle- aH e- pJ Sie~ l3H v_ , (3.6) 

where 

x = e~ 4p . (3.7) 

This corresponds to taking the limits K — * 0, j, M — j — > oo in the matrix 
elements of (I2.19|) . while keeping a = jK and (3 = K(M — j) fixed. 
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The advantage of these expressions is that they are finite-lattice partition 
functions, but with T J replaced by the rather simpler matrix e _aW . For 
definiteness, we have taken the i in (|2.19p to be 1. 

We have introduced the factor e~ p ^ into W because the elements of the 
diagonal matrix J take the values 0, 4, 8, . . . , 4[L/2], where [x] is the integer 
part of x. Hence W(a,[3,x) is a polynomial in x of degree [L/2]. This nat- 
urally manifests itself in the following working and provides a useful check 
against errors. 

If p — > oo, so x — > 0, then 



hence 
Also, 



~ pJ Si -> u + u\ - u_ul = v + v ] _ + v_v\ , (3.8) 



W(a,p,0) = Z + (a)Z_((3) . (3.9) 



W(a,0,x) = Z + (a) , 

W(Q,(3,x) = Z_((3) . (3.10) 

When a is large and positive, we can replace e _aW in the above definitions 
by e _aA+, i/'+V ; + or by e -aA -?/;_ -(/;!, according to the sub-space in which TL is 
acting. One can then verify that 

M- lim - ' T(a - ft1 ' - . (3.11) 

a,/3,L-oo [Z+(2a)Z_(2/?)] 1 /2 V ' 

We show in the following sections that we can use the techniques of 
Kaufman[2j and Yang[5j to express the Z,W expressions as square roots of 
L by L determinants, and hence evaluate M by using Szego's theorem on 
Toeplitz forms. [H [4] We further show that the determinants themselves can 
be expressed as squares of m by m determinants, where L — 2 < 2m < L. 



4 The Clifford algebra 

General remarks 

Set 

Dj = iCjSj = -iS,C 3 , (4.1) 

Tj = C1C2 ■ ■ ■ Cj-iDj , TL + j = C1C2 ■ ■ ■ Cj-iSj (4.2) 

for j = 1, . . . ,L. 

These F\, . . . , T 2 l anti-commute: 

T l T j + T j T i = 25 hJ l (4.3) 

for 1 < i,j < 2L, so form a Clifford algebra. If P is any orthogonal (not 
necessarily real) 2L by 2L matrix and 

2L 

H = zZ P hi T i » ( 4 - 4 ) 
3=1 
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then it is also true that 



r*r* + r*r* = 2S itj i (4.5) 

If V is a 2 L -dimensional matrix such that 

2L 

VTiV- 1 = Y, v w T i > ( 4 - 6 ) 

i=i 

for i = 1, . . . 2L, then we say that the 2L by 2L matrix y with elements Vij 
is the representative of V. If V is real and symmetric, then V is orthogonal 
and hermitian. 

Such matrices form a group C[10] If V,V' have the property (14. 61) . with 
representatives V,V, then so does VV', and its representative is VV. Also, 
if v = V', then 

y = c y' , 

where c is a scalar factor. This means that a representative matrix determines 
its parent to within a scalar factor. 
If Ti,T* are related by flQD, then 

VTIV' 1 = ^1^* , (4-7) 

3 

where the v*j are the elements of the matrix PVP T . Thus mapping Tj to T* 
takes V to PVP T . Similary, if a matrix 2 L -dimensional satisfies 

2L 

HTi - TiH = (H, Ti) = Y, hi T 3 , ( 4 - 8 ) 

i=i 

where the hij are elements of a matrix i?, then the same mapping takes H 
to PHP T . We call If the H-representative of ii\ 

Any representative V is orthogonal and any if -representative H is anti- 
symmetric. If V (H) is hermitian, then (because the Tj are hermitian), then 
V (H) is also hermitian. 

Exponentials of H, H 

If H is hermitian, then it must be pure imaginary and anti-symmetric. It is 
diagonalizable by a unitary matrix. Further, its eigenvalues are real and occur 
in pairs A and —A, the corresponding eigenvectors being complex conjugates 
of one another. Adding and subtracting these eigenvectors, it follows that 
there is a real orthogonal transformation that takes Ti to the form 




where D is a diagonal matrix. 

By making such a transformation and then explicitly considering the two- 
by-two constituent blocks of T>, we find that exp(— aTC) 6 Q, so 

representative of exp(— dH) = exp(-aTt) . (4-10) 
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Partition functions 

If V is hermitian, we can similarly find an orthogonal transformation (14. 4D that 
takes V to block diagonal form, each block being hermitian and orthogonal, 
so of the form 

(cosh2aj isinh2aj \ 
—l sinh 2a j cosh 2a j J 



(4.11) 



where j = 1, . . . , L. If we arrange the T*j so that this block is in rows and 
columns j, L + j, it follows that 



V = cexpL^^n+A > ( 4 - 12 ) 

where c is some scalar factor. 

Each term T^T* L+ - commutes with all the other terms, so we have decom- 
posed V into a direct product of two-by-two matrices, the jth such matrix 
having eigenvalues e aj , e~ aj . Hence 

L 

trace V/(det V) 1/d = J^2cosh ai , (4.13) 
i=i 

where d = 2 L . On the other hand, from (|4.1ip . 

L 

det(/ + t>) = n4cosh 2 ai , (4.14) 



so we have the identity 



trace V/(detV) 1/d = [det(J + V)] 1/2 . (4.15) 

Note that the relations (|4.10p . (|4.15p are quite general, being unchanged by 
the orthogonality transformation (14. 4p . 

Representative of exp(— aH) 
From ((331), (|Q> . 

L-l 

Ho = iRT L T L+1 - i £ r i r L+J+1 , (4.16) 
i=i 

L 
3=1 

so we see that Ho, TL\ are both quadratic forms in the Ti, provided we intro- 
duce the spin-reversal operator R. We can do this, since in either sub-space 
V±, R commutes with T and 7i, having the value r = +1 in V+, and r = — 1 
in V_. 
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From now on we take Hq to be defined by (|4.16p and J by (|3.4|) . with R 
fixed to be either r = +1 or r = — 1. The minimum eigenvalue of J7 is still 
zero, but is now unique, the corresponding eigenvector being either v+ or v_. 

To calculate the matrix elements Z, W using (|4.15p . we shall need to write 
each as the trace of an operator. We can do this by introducing the operator 
exp(— 7i7). In the limit 7 — > +00 we have 



v r vi 



It follows that 



(4.17) 



[Wo, Ti] = 2iTl + 



-i+i 



[Wo,r LH 



-2ir 



i-1 



[Wi.r.] = -2ir L , 



[Wi,r L 



2ir, ; 



(4.18) 
-rTi+i and 



for i = 1, . . . , L, provided that on the RHS we take I^l+i 
r = -rT L . 

We shall need some sparse L by L matrices. Let 0, I be the zero and 
identity L by L matrices, respectively, and let A,B be the one-off and two- 
off-diagonal L by L matrices 



( 1 • • ■ • \ 
• 1 



A 



\ 



-r 



1 
1 
• • 



( 



B 



10- 
10 





-r 
-r 



Note that 



AA T = A T A = I , B = A 2 
Then from (I4.18h . the -^-representatives are 



1 

• • 



(4.19) 



Ho 



2i 




A 1 



-A 




(4.20) 



W1 



2i 



/ 
-/ 



(4.21) 



Diagonalization of H 

From (13. 2p . the //-representative of TL is TLq + k'TCi. 

We calculate this matrix. First we define a 2L by2L matrix M (not to be 
confused with the M of section 2): 



M 




(4.22) 
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(writing 1 for 1/ and 1 for I). Then 



M^HoM = 2 I 1 J J , (4.23) 



M-^M= . (4.24) 



^ + ^ T / - B 
/ — £> T -A- A 1 , 

Thus M^HqM is diagonal. 
For j = 1, . . . ,L, define 

% = vr(2j - 1)/L if r = + , 

27rj/L if r = - , (4.25) 

set 

z- ~ e 1 ^ 
^ — e , 

and let P, Q be the L by L matrices with entries 

= 4"V^ , Qij = iz}' 2 /VI . (4.26) 

These matrices are unitary, but not orthogonal. Set 

P 
Q 

Then 



V = ( P ° ) • (4.27) 



U x = V~ X M- X U\MV = 2 



H = V^M^HqMV = 2 ^ ^ J j , (4.28) 

^ S -C ) ' 

where C, 5 are diagonal L by L matrices with diagonal entries 

Cjj = cos 9 j , Sjj = sm9j . (4.29) 

Thus H = Hq + fc'7l!i can be re-arranged as a matrix consisting of L two- 
by-two diagonal blocks. It is then straightforward to calculate its exponential, 
giving 

exp(-aH) = {^i ^ , (4.30) 

where U = £/ r (a), V = V r (a), T = T r (a) are diagonal matrices with diagonal 
entries (for j = 1, . . . , L) 

~ , , , > 1 COS $-1 . , . , 

c/jj = Uj{r, a) = cosh(2aAj) H sinh(2aAj) , 

X j 

V j3 = Vj (r,a) = - fc/sin ^ S ; nh(2aAj) > (4-31) 

A? 

ATI / \ 1 / \ \ ^ COS ^ -1 , » 

Tjj = tj(r, a) = cosh(2aAj) smh(2aAj) 



and 



Xj = (l-2k' cos 6j + k' 2 y /2 . (4.32) 
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Calculation of Z 

We want to calculate Z + (a) and Z_((3) from (13, 5h . First consider Z+(a), so 
take r = + in the above equations. If we are working in the sub-space V+, we 
can use (|4.17p to write (|3.5p as 

Z + (a) = trace e^e -0 * (4.33) 

in the limit 7 — > +00. Since Hq and Hi are traceless, 

det e- pn ° = det e~ aW = 1 (4.34) 

and from (|3.4p . dete~ 7 ^ = e~ Ldj , J = 7Yo- From (|4.15p we therefore have 

e 2L7 Z + (a) 2 = det(/ + e-^°e- a ^) = det(e^° + e~ a ^) . (4.35) 

Using the similarity transformation of (|4.28[> . we can replace Hq, H in (14.350 
by H , H, giving 

e^ +(Q f = detf e " 7 { + !? / .) . (4.36) 

In the limit of 7 large we expect the e _27 I, V , T blocks in this determinant 
to become relatively negligible. The factors involving 7 then cancel, leaving 

L 

Z+(a) 2 = det U = det U+(a) = Uj(+, a) . (4.37) 

j'=i 

Similarly, 

L 

Z_(a) 2 = detl7_(a) = n u i(->«) • ( 4 - 38 ) 

3=1 

Calculation of W 

The function W(a,f3,x) is defined by (|3.6p . We note that the operator Si 
takes a vector in V+ to one in V_, and vice- versa. In particular, 

Siv r = v_ r 

so 

lim e~ 7j Si = v_ r vl . (4.39) 

7— >+oo 



We can use this to write (|3.6[) as a trace: 

W(a,/3,x) = trace e-^Sie-^e-^Sie^ , (4.40) 

in the limit 7 — ► +00. We take the factors e -7 ^, e - ^ to be acting in V_, the 
two other exponential factors to be acting in V. . 

The reason we are able to make further progress is that Si also belongs to 
the group Q. In fact Si = Tl+i- This is linear in the Tj, unlike Ho and Hi, 
which from (I4.16|) are quadratic. Even so, 

SlI^ST 1 = ejT,- , (4.41) 
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where = +1 3 else e,- = —1. The representative of S\ is therefore 



5! 



-/ 
—E 



(4.42) 



where E is the diagonal matrix 



( -I • • • • \ 
1 



E 



\ 



1 

• • 1/ 



(4.43) 



Remembering that J = 7io, the representative of the matrix whose trace 
is to be evaluated in (|4,40p is therefore the 2L by 2L matrix 



W(a,(3,x) = e- 7Wo 5ie- a7 V pWo 5ie- 



PH 



and from (|4.15p 



i 2L ^ +p) W(a,f3,x) 



det 



I + W(a,(3,x) 



(4.44) 



(4.45) 



We now use the similarity transformation of (|4.28[) . We have to be care- 
ful because A,M,P,Q,V depend on r. We write them more explicitly as 
A r ,M r ,P r ,Q r ,V r - Similarly we write the 8j, Zj of (|4.25[) - (|4.27p as 9 r j, z r j. 
For the middle two exponential factors we take r = +, for the others we take 
r = — . The result is 



yy = VZ 1 MI X WMJP_ = e^™ S_e~ a ' fL e- p ' fio S+e~ m , 



(4.46) 



where 



S r = V- l M~ 1 SiM_ r V_ 



This threatens to become messy, but we find some remarkable simplifica- 
tions. Firstly, 

M r 5i = 5iM_ r , 
so M~ l SiM^ r = Si and (fOTj) becomes 

S r = V~ x S{P_ r . 
The second surprise is that 

Qr E Q — r — P r P—r — Ey , 

where F r is the L by L matrix with entries 



satisfying 



7 (r) = 

F_F+ = I . 



(4.47) 
ilifica- 

(4.48) 
(4.49) 

(4.50) 

(4.51) 
(4.52) 
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Hence S r has the comparitively simple block-diagonal form 

S r = 

Let 



F r 
F r 



Y = S_e- aH e- pHo S+e-W 
then from (ji~45j) and (14461) 



det (I + e 



-7W0 



Y 



det (< 



,7Wo 



(4.53) 
(4.54) 
(4.55) 



The RHS has the same stucture as that of (I4.36j) . so taking the limit 7 — ► +00, 
we obtain, similarly to (|4.3T|) . 



e 2Lp W(a,(3,x) 2 = detY u 



(4.56) 



where Y n is the top-left L by L block of Y. Using (p~28|) and (QUI) , it follows 
that 



W(a,/3,x) 



det 



F_ U + (a) F + U_ ((3) +xF_V + (a) F + V_ ((3) 



Set, for r = + or — , 



XJa) 



U r (a)] 1 V r (a) 



(4.57) 



(4.58) 



Then from (|4.37p and (|4.38p . remembering that the U, V matrices are diago- 
nal, so commute, we can write (I4.57P as 



W(a,(3,x) 



Z+(ayZ_(j3Y det I + xF+ l X+(a) F + X_(0) 



(4.59) 



To summarize thus far, the equations f|4.3T[> . (|4.38[> . (|4.59|> give Z + (a), 
Z-(j3), W(a,(3, x) in terms of L by L determinants, and these formulae are 
exact for finite a, (3, L. The spontaneous magnetization M is then given by 
(|3.1ip . taking the limit a, (3,L — > +00. 



5 Calculation of J\A 

So far we have parallelled the method of Yang. |5j] We now explicitly calculate 
M from the above equations, but we use a method more like that of Montroll, 
Potts and Ward.g] 

First set x = 1 and take the limit a,/3 -> +00 in (p~37j) . (HT381) . (p~59|) . 
using (|4.3ip . We obtain 



where 



M 2 = i+i_det I + F+ l X + (oo) F + X_(oo) 

1/2 



n 

3=1 



Ay } j ~\~ 1 — COS 



2A 



(5.1) 
(5.2) 



'■■j 
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and X r (oo) is the diagonal matrix with diagonal entries 

—k! sin 9 r 



From (|4"30|) . F + = -P+ 1 P_, so if we set 



Ay^' ~t~ 1 — k COS 



X r = P r X r (oo)p- x 



then (|5.ip can be written 

M 2 = £+£_det 
We could write this determinant as 



/ + X + X_ 



(5.3) 



(5.4) 



(5.5) 



detX + dot 



XT 1 + X_ 



From the definition (|4.26p of P, the matrices X+, X_ are anti-cyclic and cyclic, 
respectively, i.e. 

[X r ]i j = rji-j = -r T)L+i-j , 

where rj is defined by this equation. The matrix X+ 1 + X_ is therefore 
Toeplitz, its elements i,j depending only on i — j, and one might hope to 
use Szego's theorem^ |8] to evaluate the determinant for large L. 

Unfortunately this naive approach does not work because for large L the 
determinant is not dominated by its near-diagonal elements, in particular 
elements such as (L, 1), (2, L— 3) do not tend to zero. However, we can rescue 
this idea by making a particular bilinear transformation of X + , X_ . 

We set 

C r = {I -iX r y l (I + iX r ) = P r c r P~ l , (5.6) 
where c r is a diagonal matrix with diagonal elements 



( c r)j,j 



X r j + 1 — k' exp(i0 r j) 
X r j + 1 — k' exp(— iO r j) 

Solving for X r and substituting into (|5,5p . we obtain 
M 2 = eU-det[(C + + C_)/2] , 



(5.7) 



(5- 



where 



C = £ r det(I-iX r 



n 



[A rJ + l-fc / exp(-ig r , i )] 2 
^ 1 2 X r> j (X r j + 1 — k' cos 6 r j) 



1/2 



(5- 



From (ET32j) . 



[1 - k' exp(i9 r>j )] [1 - fc / exp(-i6> r . !i )] 
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from which it follows that 



e = n 



1 — k' exp(— i6 r ,j) 
1 — k' exp(i6 r j) 



1/4 



(5.10) 



The 9 r j are either ir or 2ir, or occur in pairs 9 r j,2n — 9 r j. It follows that 



Also 



simplifies to 



1 — k' exp(—i6 r j) 
1 — k' ex.p(i6 r j) 

From (|4.26p and (|5.6p . the elements of C r are 

1 



.1/2 



(5.11) 



(5.12) 



(a 



L 

E 



^ / j \-t Jm,m^ r ,m ) 
m=l 



(5.13) 



where z r 



The summand depends on ?Ti via r m and is an analytic 



function of 9 r m on the real axis, periodic of period 2ir. The 9 r m are distributed 



uniformly throughout this period. If i,j are held finite and L 
that 



where 



f(0) 



2tt Jo 
1 — k' exp(— i 



1 



1/2 



oo, it follows 
(5.14) 

(5.15) 



k' exp(i^) 

This is actually the function e 1<5 * of OnsagerfTJ eq. 89] and Montroll et aZ[H 
eq. 42] In their notation z\ = tanhff , z| = tanh^T* and our k' is k' = z^/zi. 
We have replaced u by 9 and taken the hamiltonian limit Zi, z\ — > 0. 

Note that the limit (|5.14p is the same for r = +1 and r = — 1. Also, for 
large but finite L, the corrections to (|5,14p vanish exponentially with L. The 
near-diagonal elements of C+ and C_ therefore become equal when L is large. 

Since z^ m = —r, incrementing i or j in (I5.13P by L negates (C r )i,j for 
r = +, and leaves it unchanged for r = — . This means that the near-top-right 
and near-bottom-left elements of C+ and C_ also tend to finite limits as L 
becomes large, but are equal and opposite. Their sum therefore approaches 
exponentially to zero . 

Hence when L is large we expect (|5.8p to become 



detC 



(5.16) 



where C is the L by L matrix with elements given by (|5.14p . These 

elements tend exponentially to zero as \i — j\ becomes large. This is a special 
case of the Toeplitz matrix discussed in [3] and we can use the general result 
(68) therein: 



lim G~ L det C 



exp 




(5.17) 
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where 

G = exp(^jf^log/(0)dl?) , 

and 

oo 

io g /(0) = ^ md ■ 

n=—oo 

Since log f{9) is an odd periodic function of 6, we have G = 1. Also, taking 
the logarithm of (|5. 15[) and Laurent expanding, we readily find 

2 Kn = k' n /n if n > , 2n n = -k'~ n /n if n < . (5.18) 

It follows that detC = (1 - fc' 2 ) 1 / 4 , and hence 

M = (1-k' 2 ) 1 / 8 . (5.19) 

This is of course Onsager's famous result. [H O [5]. 

6 Connection with the superintegrable chi- 
ral Potts model 

The superintegrable case of the iV-state chiral Potts model has some properties 
that closely resemble the Ising model. In particular, there is a spin-shift 
operator R that is the natural generalization of (|2.14p and divides the vector 
space into N sub-spaces Vq, labelled by Q = 0, 1, . . . N±. Within each Vq, if 
one imposes the fixed-spin boundary conditions of Fig. 1, the transfer matrix 
T generates a yet smaller sub-space of dimension 2 m , where 



m = rriQ 



(N -1)L-Q 



N 



(6.1) 



Here [x] means the integer part of the real number x. 

Within Vq, T is a direct product of m two-by-two matrices, and there 
are similarity transformations that reduce the associated hamiltonian to the 
direct sum[TTl eq.2.20] 

in 

H = n Q - N [(1 - k' cos Oj)Sj- k! sin 9jCj] , 
i=i 

where 

fj, = fj,Q = 2k' Q + (1 + k')(mN - NL + L) (6.2) 
and 0\,...Q m are defined by 

cosOj = (1 + Wj)/(1 - wj) , < 0i < 7T , (6.3) 

the u>i, . . . , w m being the zeros of 

N-l 

P(z N ) = Z -Q J L+Q)n [(z N - \)/{z - co n )] L , (6.4) 

n=0 

which is a polynomial in w = z N of degree m. 
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The partition functions Z 

One can explicitly calculate [12] the partition functions that generalize (j4.37j) . 
(PSD : 

Z Q (a) = e-^ Ul {a)---u m (a) , (6.5) 

the function Uj(a) again being defined by (|4.3ip . (|4.32p . 
For N = 2 the model reduces to the Ising model with 

r = 1 — 2Q , m(r) = tuq , (6-6) 

so Q = corresponds to R = +1, while Q = 1 corresponds to R = — 1. The 
result (|6.5p of course agrees with (|4,37p . (|4.38p . but has a slightly different 
form. The partition function is no longer squared on the LHS and instead of 
there being L variables 6j, there are only m, which lie between (L — 2)/2 and 
L/2. They are the same as the Oj of the previous sections, but j takes only 
the values 1 to m and < Oj < tt. 

These differences are easily explained. For every Oj with j = l,...,m 
there is a Qy = 2n — Oj, where j' = L + 1 — Q — j and m < j' < L. They 
have the same value of Uj(a) and each occurs in (|4.37p . (|4.38p . while only one 
occurs in (16. 5p . This accounts for the absence of the square in (16. 5p and the 
presence of the factor e~^ a accounts for the exceptional cases Oj = tt or 2tt, 
which are included in the product in (|4.37|) . (14.381) . 

There are a total of four cases to consider. 

1) r = +1, Q = 0, L = even: there are no exceptional cases and m = L/2, 
/x = 0. 

2) r = +1, Q = 0, L = odd: then m = (L — l)/2,/x = -1 - k'. There 
is an exceptional case at j = (L + l)/2 = m + 1, where Oj = tt, Uj(a) = 
exp[2(l + k')a] = e" 2 ^. 

3) r = -1, Q = 1, L = even: then m = (L — 2)/2,/i = —2. There 
are two exceptional cases at j = L/2 and L, with 8j = tt and 2ir, Uj(a) = 
exp[2(l + k')a] and exp[2(l — k')a], so u L / 2 (a) ul(o) = e~ 2 ^ a . 

4) r = -1, Q = 1, L = odd: then m = {L — l)/2,/x = k' — 1. There is an 
exceptional case at j = L, with 6j = 2tt, Uj(a) = exp[2(l — k')a] = e~ 2>Ma . 

In each case the total contribution of the exceptional cases to the RHS of 
(1437]) or (f4~38|) is e~ 2/iC \ which accounts for the factor e~^ a in (163]) . 

The partition function W 

The equation (|6.5p simplifies (14.37P and (|4.38p . expressing the Zq{q) as a 
product rather than the square root of a product, and reducing the number 
of factors from L to m or ra', where 

m = m(+) = [L/2] , m! = m(-) = [(L - l)/2] . (6.7) 

Can we similarly reduce the expression (|4.57p for W(a,{3, x)? 

The answer is yes, apart from simple factors that are independent of x and 
easily calculated. We show that the matrix 

G = U + {a)F + U_{l3) + xV + {a)F + V_{f3) (6.8) 
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can be reduced to block lower-triangular form and its determinant expressed 
as simple factors times the square of an m' by m! determinant. 

For brevity, in the following two sub-sections, unless indicated otherwise 
we write 



9j - 9+j , 6j - 9-j , 



u 



U+(a)j,j = u +:j {a) , Vj = V+(a)jj = v + j(a) , 
= U-(0) j j = u- d (P) , v'j = V_(f3)jj = v—j({3) . 



We use the definition (IPT1) of F± and the properties (|3~52|) . fTT2l) . We 
also implicitly use 

UL+l-j = Uj , VL+l-j = -Vj , u' L _j = Uj , v' L __j = -v'j . (6.9) 

We do not in this section need the definitions (I4.31|) . 

The case L even 

In this case note that the U+(a), V+(a) on the left of (j6.8jl depend on the 
9j, which satisfy 9i + 9i+i-i = 2ir, so it is natural to combine rows i and 
L + l-i. The U_{{3), V_(fl) on the right involve 9'j, satisfying 0'j + 9' L _j = 2vr, 
so it is natural to combine columns j and L — j. Also, columns L/2 and L 
both correspond to exceptional values of 9'j, where sin^- = 0. The elements 
v 'l/2i v 'l °f vanish for these two columns. From (|6,7I ). 

m = L/2 , m = (L-2)/2 = m-l . (6.10) 
We perform the following equivalence transformations sequentially on G: 



1) Gij - 


G^ + Gi+i-ij, 


1 < i < m, 


2) Gij — * 


Gij ~\~ Gi L—j, 


m < j < L 


3) G hL - 


— UjJGi^ m l 'u m , 




4) - 


Gij VijGi^ m /u m , 


1 < j < 771 


5) Gjj — i 


Gij 2UjGj rn /'u rn , 


m < j < L 


6) Gij —s 


■ Gy-(l + ^>5G iiL K, 


m < j < L 



(6.11) 



where = cos(0j). 

The first transformation corresponds to pre-multiplying G by some matrix, 
all the others to post-multiplying it. Steps (2) to (6) are for all values of i, 
i.e. 1 < i < L. At every step we are merely incrementing rows or columns 
by linear combinations of other rows or columns, so the determinant of G is 
unchanged. 

The elements F^j of F + satisfy 

Fij + F L+1 ^ L ^ = F i)L + F L+1 ^ L = -2/L (6.12) 
for all It follows that the final matrix G is such that 

dj = for 1 < i < m, m<j<L . (6.13) 
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Thus G is now block lower-triangular: 



where the Y, T, Z are m by m matrices. We find 

-2i (sin OjUiUj + x sin QiV{v\ 
L(cos #j — cos 6j) 



for 1 < i < m, 1 < j < m. 

Note that the equivalence transformations are independent of x, so Y, T, Z 
are all linear in x, their coefficients being the transforms of the corresponding 
coefficients of G. 

We define the m by m matrix Z by (|6.14j) , but with the rows and columns 
re- arranged so that, for 1 < i,j < m, 

(Zh = (G)L + i-i,f , (6-16) 

where f = L if j = m, else f = L — j. Let D, D' be the diagonal matrices 
with elements 

D jtj = sm6i , D' j:j = sinO'i (6.17) 
for j = 1, . . . , m, except that D' m m = —iu' L /u' m . Then we find that 

Z = —D~ X YB' . (6.18) 

Hence 

det G = det(-D')(det Y) 2 / det D (6.19) 

and we see that det G is basically the square of the mbyra determinant det Y. 

We can do better yet and write it in terms of the square of an m' = m — 1 
determinant. Expand the matrices G, Y in powers of x: 

G = g + xgi , y = y + xyi , (6.20) 

so go,gi,yo>yo ar e ^-independent matrices and go = U + (a) F + U_((3). From 
(|4.52p . its inverse is 

go 1 = U_(P)- 1 F_U + (a)- 1 . (6.21) 

We can follow the above six steps (|6.1ip and evaluate gp 1 after the six equiv- 
alence transformations. It also becomes block lower-triangular, as in (|6.14j) . 
In particular its top- left block is yo 1 ) with elements 



yo 1 



-2isin# 



/.(( OS 0' - cos Q^u'jUj 



yo 1 . = -VKn«j 
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for 1 < i < m, 1 < j < m. This matrix is necessarily the inverse of yo and 
from (I6TT9]) 

detg^G = [detyo 1 ^ 2 . (6.23) 

We calculate y^ 1 !". From (I6.15h . the last column o/yi vanishes. Hence, 
using (|fT20jh 

y -!y = I + X f) , (6.24) 

where T> = y$ 1 yi and again the last column of T> vanishes. The determinant 
of the RHS is therefore the same as that of the RHS truncated to its first 
m! = m — 1 rows and columns. 

We can therefore ignore the last column of Y and the last row of y^" 1 when 
calculating the LHS, and contract Y to the non-square m by m' matrix with 
elements given by the first of the equations (|6.15p . y^" 1 to the m! by m matrix 
with elements given by the first of the equations ()6.22[) . 

We define the m by m! matrix A + and the m' by m matrix A_ by 

/a \ 2 sin6*j , . 2 sin^ . . 

{A+kj ~ L(cos^-cos^) ' ( kj ~ L(cos^-cos^) ' [ ' 

Then from (pTKj) . ^22\i . 

yi = -iV+(a)A+V_(P) , y^ 1 = -iU-{p)- x A_U + {a)- 1 , (6.26) 

where we now take the U r , V r ,X r to be the diagonal matrices defined by (|4.3ip . 
(|4.58p . but truncated to the first m(r) terms. 
Hence 

v = -u^r l vu_{0) , 

where 

V = -A_X + {a)A + X_(fi) . (6.27) 
From (I6.23p . we obtain 

detG = det (u+ (a) F + U_(/3)) [det(I m , + xV)} 2 , (6.28) 

where I m i is the identity m! by m! matrix. 

The relations (|4.57p . ()4.59p therefore become 

W(a,fi,x) = Z+(a)Z_(p)det[I mf - xA_X + (a)A + X_(f3)} . (6.29) 

We have eliminated the square on the LHS and used only matrices of dimen- 
sions m or m', involving only 0\, . . . , 6 m , 6[, . . . , 0' m ,. 

The case r = +, L odd 

Now (|6.7p gives m = vnl = (L — 1)/2, so the matrix G has one exceptional row 
at i = m + 1 and one exceptional column at j = L. The steps corresponding 
to (I6TTT|) are 

1) Gij -> Gij + G L +i-i,j, 1 < i < m, 1 < j < L 

2) G^ G^ - 2uiG m+ ij/u m+1 , l<i<m, l<j<L 
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3) Gij - 

4) Gij 
5) Gij - 



Gi.4-2u' j G i , L /v? L , 

4 



G 



UjG itL /u' L 



m < j < L 
m < j < L 
1 < j < m 



(6.30) 



The first two transformations correspond to pre-multiplying G r by some ma- 
trix, all the others to post-multiplying it. Steps (3) to (6) are for all values of 
1 < i < L 

Using (|6.12p . we find that the transformed G has the structure 



G 



( Y \ 

•• -u m+ iu' L 

\T z •• j 



(6.31) 



where Y, T, Z are all m by m' = m blocks and the middle row and last column 
have width one. Hence 



we also find 



Yi 



'j 



detG = u m+ iu L (det Y) (det Z) 



-4i cos 2 (OiWluiu'j tan(0j/2) + xv^ tan(0;/2)] 

cos B'i 



(6.32) 



/.((•OS !., , 

and again if we rearrange Z so that Z%j = Gl+i-i^l-] for 1 < i, j < m, then 

Z = -D~ X YD' , (6.33) 
with D, D' given by (|6.17p for all 1 < j < m. Further, 

-2isin6>' 



yo 1 )^- 



L(cos d\ — cos Qj^u'jUj 



(6.34) 



From (|6.15p and (|6.32|) . yo is different for the two cases (L even and L 
odd). However, the wanted elements of yi and of y^ 1 are the same, so again 
W(a,(3,x) is given by (j6T29]) . with A+, A_ given by (f6T25|) . 

From ()4.3ip . (|4.58p . the diagonal matrices X± have diagonal elements 



[*+(«)] 



k' sin 6j sinh(2aAj 



jj 



Xj cosh(2aAj ) + (1 — k' cos 6j) sinh(2aAj) 
k' sin^.sinh(2/3A^) 



(6.35) 



jj 



X'j cosh(2/3A^) + (1 - k' cos d'j) sinh(2/3A^) ' 
Xj being defined by (|4.32p and A^- by the same equation with 9j replaced by 
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7 W in terms of orthogonal matrices B +1 B_ 

In this section we take the result (I6.29P and write it in terms, not of the 
matrices A + , A_ , but of related orthogonal matrices B + , B_ . 

In both the L even and odd cases, the matrices Vq , y^ above have a very 
similar structure, differing only by diagonal equivalence transformations. We 
can therefore construct orthogonal matrices B + , B_ that differ from A + , A_ 
only by such diagonal equivalence transformations. 

Let E, E' be diagonal matrices with diagonal elements 

Ejj = sin 9j , E'aa = sin Q'- if L = even , 

Ejj =tan(0j/2) , E'^ = tan(^/2) if L = odd. (7.1) 

Then whether L is even or odd, we find from (16. 15ft . (|6.17p . (|6.22p . (|6.25jl . 
(16321) . (lOill that 

yo 1 =iD'A\D- x , y = -iE- 1 A + E' , A T _ = -D- X A+D' , (7.2) 

yo 1 being the left-inverse of yo- 

Hence we can define matrices B + , B_ : 

B + = {DE)- 1 / 2 A + {D'E') l l 2 , B_ = {E' / D') 1 ' 2 A_{D / E) l > 2 (7.3) 

such that 

B_B T = BlB + = I m , , B_ = -Bl , (7.4) 

so, even though they are not necessarily square, they are orthogonal in this 
sense. Their elements are 

2 sin 0'- 

(B +)lJ = -(fl_) i<= 7^73^ if ^= even, 
/ p % , R v 4cos(^/2)sin(0;./2) 

(5+)ii= -(B_)i<= pTTTTTi ZTT^n lf L = odd ■ ( 7 - 5 ) 



L(cos 9i — cos 0j) 



Set 



E + = E , ^_ = E 1 ' -1 . (7.6) 
then we can write the result (|6.29[) as 
W(a,P,x) = Z+(a)Z_(f3)det[I m > - xX + (a) E + B + X_(f3) E^B_] . (7.7) 

We can slightly generalize this result. Write W(a,f3,x) as W + (a, (3,x), 
and define a new function W_(a,/3, x) to be given by (|3.6[) with v+,v_ inter- 
changed. Thus (|3.5p . (|3.6p become 

Z P (a) = w|e aW x; P , W r (a, P, x) = vle aH e pJ S^v^ . (7.8) 

Then we can generalize (|7.7p to 

W r (a,p,x) = Z r (a)Z_ r (P)det[I-xX r (a) E r B r X_ r (P)E_ r B_ r ] , (7.9) 

for r = ±, where I is the identity matrix of dimension m(r). For r = +, 
(|7.9I) is (I7.7p . For r = — it can be deduced from it by taking the hermitian 
conjugate of the definition (|7.8p of W_ and interchanging a with p. We also 
need the general identity 

det(J + AB) = det(7 + BA) , 

which is true for all matrices A, B even when they are non-square, 
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8 Conclusions 



We have considered the zero-field Ising model on a cylindrical lattice with 
fixed spin boundary condtions on the top and bottom rows, and have replaced 
the transfer matrix product by the exponential of the associated hamiltonian. 
This leaves the eigenvectors and spontaneous magnetization unchanged. 

We have then used the Clifford algebra technique of Kaufman [2] to evaluate 
as L by L determinants the partition function Z, and the partition function 
W with a single-spin operator Si included. In section 5 we show how to 
use Szego's theorem on Toeplitz matrices to then obtain the spontaneous 
magnetization. 

Much of this merely parallels Yang's derivation [5], but in section 6 we go 
on to reduce Z and W to determinants of size m or m' where L — 2< 2m < L. 

In this form we can compare the Ising model results for Z with those of 
the iV-state sup erintegr able chiral Potts model, taking the hamiltonian limit 
therein. [TTlll2j The Ising model is the superintegrable chiral Potts model with 
N = 2, and the results of course agree. 

All this leads up to the question whether we can generalize the Ising result 
(|6.29|) or (|7.9p for W to the Y-state superintegrable chiral Potts model, for 
general N. This should open up the possibility of an algebraic derivation of 
the spontaneous magnetization of that model. We have numerical evidence 
that strongly suggests the answer to the question is yes, and that one can 
make a fairly immediate generalization of the form of the result given in (17. 9ft 
in terms of the orthogonal matrices B + ,B_. We shall present the conjecture 
in an accompanying paper. [12] 
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